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We study double transverse spin asymmetries in the Drell-Yan process at measured transverse 
momentum of the lepton pair. Contrary to what a collinear factorization approach would suggest, a 
nonzero double transverse spin asymmetry in the laboratory frame a priori does not imply nonzero 
transversity. TMD effects, such as the double Sivers effect, in principle form a background. Using the 
current knowledge of the relevant TMDs we estimate their contribution in the laboratory frame for 
Drell-Yan and W production at RHIC and point out a cross-check asymmetry measurement to bound 
the TMD contributions. We also comment on the transverse momentum integrated asymmetries 
that only receive power suppressed background contributions. 
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I I. INTRODUCTION 

Transversity -the distribution of transversely polarized quarks inside a transversely polarized hadron- was first dis- 
cussed by Ralston and Soper who suggested its measurement in the double polarized Drell-Yan process. Although 
this suggestion was made more than 30 years ago, this demanding double transverse spin asymmetry measurement 
has not yet been performed. At present it is in the future physics program of BNL's Relativistic Heavy Ion Collider 
(RHIC) 1 and is also considered at GSI-FAIR, J-PARC and NICA. Currently RHIC is the only accelerator where 
polarized proton-proton collisions can be performed. Therefore, in this article we will focus on RHIC. 

Ralston and Soper considered the double transverse spin asymmetry (Axx) integrated over the transverse momen- 
tum qx of the lepton pair (later also reconsidered in [SriEj]), and at measured qx, in particular at qx = 0. Both cases 
£f~) | will receive nonzero contributions from transversity. The asymmetry as a function of qx has been studied in Ref. @ 
. in a collinear Collins- Soper- Sterman resummation approach, showing it to be maximally of order 5% and fairly fiat 
in qx up to a few GeV. At measured qx there will however be background contributions from transverse momentum 
dependence of partons, that have not yet been considered. We will study these contributions to the double transverse 
spin asymmetries at measured transverse momentum in both the Drell-Yan process and in VF-boson production, where 
in the latter case one expects zero contribution from transversity within the Standard Model @, ■ 

At RHIC single transverse &phi asymmetries (An) will be studied in W production as well, with the goal of measuring 
the sign of the Sivers effect [U[l(J. This effect refers to the fact that the transverse momentum distribution of quarks 
inside a transversely polarized hadron can be asymmetric w.r.t. the spin direction [Tl| . This spin-orbit coupling 
effect arises from initial and/or final state interactions and has a calculable process dependence (when factorization 
applies). The sign in Drell-Yan or W production is predicted to be opposite to the one in semi-inclusive deep inelastic 
scattering (SIDIS), the process in which the Sivers effect asymmetry was first observed [l2l - fl4j ]. The Sivers effect 
may also generate background for the transversity double transverse spin asymmetries in the Drell-Yan process at 
measured qx, through a Sivers effect in both incoming hadrons. This double Sivers effect will be investigated in this 
paper. Moreover, it can lead to a nonzero result in W production, which could be mistaken for physics beyond the 
Standard Model, for instance from the complex mixing of W bosons with a hypothetical W boson that appears in 
many extensions of the Standard Model [la ]. We will study these aspects of the double Sivers effect contribution 
quantitatively in this paper. 
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Besides the double Sivers effect we will include contributions from another transverse momentum dependent effect 
that was first discussed by Ralston and Soper [l| : it describes the distribution of longitudinally polarized quarks inside 
a transversely polarized hadron. Both effects are described by a transverse momentum dependent parton distribution 
(TMD): the Sivers effect by a TMD often denoted by f^p [16[ and the other by gyr Ujj- The latter function also 
appears in the analysis of the evolution equation of the twist-three function gx = gi + gi [18l422j and is in the literature 
sometimes referred to as one of the two "worm gear" functions [23| . We will adopt this convention and refer to grr 
as the Worm Gear (WG) function, because the other worm gear function h\ L will not be discussed here. 

The expressions for the double Sivers and WG effects for Drell-Yan have been given in Ref. |24l426| . Quantitatively, 
both effects have been studied in Ref. [25| in polarized proton- antiproton Drell-Yan at GSI-FAIR, and the double WG 
effect in W and Z boson production at RH1C in Ref. [27|. The latter study contains some mistake in its Eq. (21) that 
will be corrected here, without altering the conclusions. 

As can be seen from the expressions in Ref. (23 - [26| . one can consider a specific frame, the so-called Collins- 
Soper frame, that in principle allows one to distinguish the double transverse spin asymmetries Att(qt) arising from 
trans versity, the Sivers effect and the WG effect. Different angular dependences can be projected out allowing to single 
out a specific contribution. Transversity leads to a spin asymmetry proportional to cos 2c/) l s , where the azimuthal angle 
4> l s is measured between the spin plane and the lepton plane, whereas the other two effects lead to spin asymmetries 
independent of this lepton azimuthal angle. However, in the laboratory frame, as we will show in this paper, all 
three effects will contribute to the angular distribution cos 2(p l s . This is in contrast to what a collinear factorization 
approach would suggest, e.g., in the treatment as applied in Ref. the Sivers and WG contributions are absent 
from the start. In that approach the expression for Att(qt) in the lab frame will be only in terms of the transversity 
distribution (see also Ref. [l7]), which might lead to wrong conclusions about transversity. 

The lab frame is thus a priori not the right frame to extract the transversity distribution from the asymmetry 
Att(qt)i however, it is experimentally more 'direct' to extract spin asymmetries in the lab frame. Analyzing the data 
in the lab frame might also be more accurate, because any additional uncertainties from the transformation to the 
CS frame will avoided. Furthermore, an analysis in the CS frame and the lab frame could be cross-checked with each 
other if one knows the expected differences between the two frames. In VF-boson production with a leptonic decay, 
on the other hand, it is impossible to transform to the CS frame, because the neutrino will go unobserved rendering 
it impossible to determine the transverse momentum of the W boson. This means that the double transverse spin 
asymmetries in VF-boson production have to be studied in the lab frame, where the transverse momentum dependent 
effects form a background for the new physics studies as proposed in [l5j . We think it is therefore important to know 
quantitatively the size of the spin asymmetries in the lab frame caused by partonic transverse momentum effects in 
the Drell-Yan process and IF-boson production, which will be explored in the remainder of this paper. 



In both processes we have to deal with vector-boson production from hadron-hadron collisions, with a subsequent 
leptonic decay. The cross section for such a process has its leading contribution coming from the diagram in Fig. [T] 
This contribution can be split into a lepton and hadron tensor connected by the appropriate vector-boson propagator. 



II. DRELL-YAN CROSS SECTION IN TMD FACTORIZATION 
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FIG. 1: Leading diagram in the Drell-Yan process. 



The hadronic tensor can be expressed by 



W^ u = - 




(1) 
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in terms of the fully unintegrated quark correlator 



1 



(27T) 



d 4 £e tp < (P, S\qj (0)W(0, £)« (£)|P 5 ) 



(2) 



and the quark- vector-boson vertices V^, and V^/, where the primed vertex has the complex conjugated coupling 
strength. The gauge link U(0, £) that renders the correlator gauge invariant is not specified at this stage. Writing 
the momenta in terms of lightcone and transverse components as p — \p~ , p + , pt] , the delta function in the hadron 
tensor can be approximated by 



S 4 (p + k-q)^ 6(p+ - q+)5(k- - q-)5 2 {p T + k T - q T ) 



(3) 



which is accurate up to p~k + /Q 2 corrections, where Q 2 — q 2 is the photon momentum squared. Keeping higher 
order corrections in the kinematics only gives very small corrections, see e.g. [Hj]. This fixes q + — p + = XiP^ and 
q~ = fc~ = X2P2 and allows us to write the hadron tensor as 



1 
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3 ^ 



J d 2 p T d 2 k T S 2 (p T + k T - q T )Tr 



+ (lo2). 



(4) 



in terms of the p -integrated quark correlator 

S?(a;,pT)ii= / dp"$«(p,Pi,Si) 



p+=xP+ 

^^e^iPuSrimmoqmPi^i) 

and a similar $2, which is integrated over k + . The gauge link is process dependent, leading to a Sivers function in 
Drcll-Yan that has the opposite sign compared to the one in SIDIS, cf. for instance Ref. [9;] and references therein. 
In the double Sivers effect this sign is not relevant however. The quark correlators can be parameterized in terms of 
transverse momentum dependent distribution functions by 



$1(2:1, Pt) = -yi{x u p T )7/i + + f 1 ^(x 1 ,p T ) — + 9i T (xi,PT)— ^ 75#+ +•• 

frit 11 ^fflt 1 \U 1 1 \ e ^pal^n u _k^S2 T q , , -.k T -S 2 T , 1 

$2(2:2, k T ) = ^yi{ x 2,k T )i_ + f 1 ji(x2,k T ) — lgl T (x2,k T ) — — 7 5 ?*_ H j, 



(6) 



in which we only wrote the relevant distribution functions for the transverse spin asymmetries to leading order in 
M p /P^~ and M p /P 2 ". For further details we refer to, e.g., Ref. (24j |. 

We will define a spin flip symmetric and antisymmetric cross section by 

da s = - (da n + da n + do** + d<j u ) = —W^D w D* va L p °dP, 

\ ] s (7) 

da A ee ~ (da n - da n - do* + da u ) = ^-W% / D w D* l , <T L>" r dP, 

where Wg" A is the hadron tensor symmetrized or antisymmctrized with respect to the proton spins, D^ p is the 
vector-boson propagator, L pa is the lepton tensor 

L pa = Tr [VftVfl] (8) 

in terms of the lepton- vector-boson vertex Vf and, finally, dP is the phase space element 
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III. DISTRIBUTION FUNCTIONS 

The distribution function fl(x, kr) describes the probability of finding a quark q with lightcone momentum fraction 
x and transverse momentum with length kx = |kr|- As often done, for our phenomenological studies we will assume 
factorization between kx and x dependence and assume a Gaussian dependence on kx, i.e. 

fl{x,k T ) = ^ ) e- k - /{k - ) fl{x). (10) 

Such a Gaussian dependence has been shown to work very well [29] . We will use the value of the width, 

(k T ) = 0.25 GeV 2 , (11) 

found by L'-sj based on the Calm effect in unpolarized SIDIS. Although this value may differ from the (k T ) in Drcll- 
Yan, the deviation is not expected to matter for our purposes and to fall within the error in the estimates we will 
consider. 

The Sivers distribution function describes the correlation between the partonic transverse momentum and proton 
spin direction. The probability of finding a quark q with transverse momentum inside a transversely polarized 
proton is given by 

Ik IIS I 

P'fofcr) - ff(x, k T ) + sin(0 kT - c/> St ) 1 r| ' T| fer). (12) 

In SIDIS there are clear experimental observations of the asymmetries that would arise from the Sivers effect, offering 
strong support for the latter effect. Within that picture the current experimental data allows for a determination of 
the Sivers function for both the u and d quarks and anti-quarks. In the recent extraction obtained by [30l |. the Sivers 
function for SIDIS is parameterized as 

f#(x, k T ) = -Af q (x)ti(k T )f?(x, k T ) (13) 

with 

h'(k T ) = V^^e- k ^ M l 
N q {x)=N q x««{l-xf ^ + iy . 

Otq Pq 

The numerical values found in the extraction are M\ = 0.34 GeV 2 for the flavor independent width of the distribution 
and the numbers in table U for the parameters in the flavor dependent function that describes the x dependence. The 
current knowledge of the Sivers function at small x is limited, but the single spin asymmetry measurements at RHIC 
will certainly improve this. For the moment, we take what is known until a better determination will be available. 
Taking into account the error bars in [3(| we come to the rough estimate that the overall effect, as will be calculated 
in Section IIVI and El can be maximally enhanced by an order of magnitude. As said before, the sign of the Sivers 
function for Drell-Yan is supposed to be opposite to the one for SIDIS, however in the double Sivers effect this has 
no influence. 





u 


u 


d 


d 


Ctq 


0.73 


0.79 


1.08 


0.79 




3.46 


3.46 


3.46 


3.46 




0.35 


0.04 


-0.9 


-0.4 



TABLE I: Numerical values for the parameters in the Sivers function from [3C 



The Worm Gear distribution g^ T (x, kx) describes the longitudinal polarization of quarks with transverse momentum 
k<r, inside a transversely polarized proton. A determination of this distribution based on fits of experimental data is 
not available. Data on double transvers e sp in asymmetries Alt that receive contributions from the WG effect has 
become available only very recently [3lL |32| |. The recent measurements on 3 He indicate that gyr for the up-quark is 
not small 1321. 
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Here we will employ a model for this WG function. Both the bag model [33J] and the spectator model [34| agree quite 
well with a Gaussian approximation of the transverse momentum dependence for not too large values of the trans- 
verse momentum. We will therefore use the Gaussian Ansatz, which allows us to express the transverse momentum 
dependent distribution as 



9 q M*k T ) = -*-e-*/M™ g $\x), (15) 

in terms of its first transverse moment g^p (x), which is defined as 

g$\x) = J d 2 k T ^gl T ( X ,k T ). (16) 

For the width we will take a value in accordance with the bag model 

(4) wo = 0.71(fc|). (17) 



For the first moment, we will use a Wandzura-Wilczek (WW) type approximation [351 - 4371 ] to express it in terms of 
the known helicity distribution g\ [x] by 

9$${x)nx fdy 9 ^. (18) 

j x y 

For numerical estimations of this function the DSSV helicity distribution [38| will be used. Deviations from the WW 
approximation can be considered , but the WW distribution is in fair agreement with the bag model, the spectator 
model, the light cone constituent quark model [40| and the light cone quark-diquark model [4lj . Furthermore, a recent 
determination of target transverse spin asymmetries in SIDIS |42J is consistent with the theoretical prediction based on 

the WW type approximation of (43|. With all these ingredients the lowest Mellin moment gyp = I dxgyp {x) of the 
first transverse moment g^p (x) can be calculated. We find g^p = 0.091 and gfp' 1 ^ = —0.026. This is in excellent 
agreement with the evaluation on the lattice (at the scale 1.6 GeV): g r ^°' 1 ' ) = 0.1055(66) and gfp = —0.0235(38) 
from Ref. [44J. All this gives us confidence that the estimates below are sufficiently realistic. 



IV. SPIN ASYMMETRIES IN THE DRELL-YAN PROCESS 



In the Drell-Yan process the virtual photon produces a lepton and anti-lepton, both of which can be detected. This 
allows a full determination of all the kinematic variables, such that one can transform to the so-called Collins-Soper 
frame [45| . In that frame a correlation between the lepton transverse momentum and the proton transverse spin 
direction will come solely from the transversity distribution, which makes this process very suitable for an extraction 
of this distribution function from Att(qt)- If one just analyzes the correlation between the lepton angle and the 
proton spin direction in the lab frame, there can be a residual asymmetry coming from double Sivers and WG effects 
for the following reason. The Sivers and WG function both cause the photon transverse momentum and the proton 
spins to be correlated. When the virtual photon decays, the decay products are more inclined to move in the direction 
of the parent particle which, in turn, causes the direction of the decay products to be also correlated with the proton 
spin directions, albeit diluted. 

In order to estimate the error that one would possibly make in the extraction of the transversity distributions from 
Att{<1t) by not going to the Collins-Soper frame, we will calculate the double transverse spin asymmetries in the lab 
frame coming from the Sivers and WG functions. 

In the following analysis we will work towards an asymmetry differential in the photon's momentum squared 
Q 2 , transverse momentum length qx, and rapidity Y = ^\ogq + /q^. The other kinematic variables, which will be 
integrated over, are <fi q , which is the azimuthal angle of qr, and y, which is defined as y = l~ /q~. The final kinematic 
variable is the direction of the lepton transverse momentum <j>i in the lab frame, which in the end will be integrated 
over with particular weights to select out the different contributions to the spin asymmetries. 

Starting from Eq. (|7|. we can express the cross section as 

do = f ' qr ( qr cos(^ - <j> q ) 

dQdq T d<f>idY J V ^8(2^ sQ*{ ^J Q2 i_y _ # _ ^ 



\W^L 



(IV ■ 



(19) 
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by using the fact that effectively the photon propagator D^ u — —ig^/Q 2 and the phase space element in these lab 
frame coordinates is 



dP 



Qqr 

4(2tt) 



-dqi 



kdQdYdy 



q T COS(0; 



r) 



The vertices in the lepton and hadron tensor are, for the Drell-Yan process, given by 

Vgqi iCqCy Sqq', 



(20) 



(21) 



Furthermore, in the expression for the lepton tensor we need the lepton (I) and anti-lepton (I) momentum 4-vectors, 
which can be specified in terms of lightcone and transverse components, in the lab frame by 



V2 V 1-1/ 



I = 



1 

71 



1 - V -Y 

e 



(1 -y)l 2 T + yl 2 , 



{l-y)l 2 T + yl 2 r , 



2 p Y 

y ~T e 



(1 - y)l 2 T + yl\ 

J_ / V 72 p Y 
V2 V l-V l T e 

1 - y)l 2 T + yl\ 



where 



and 



q T + l T — 2lrqr cos(0i 



l T = q T ycos((f>i - 4> q ) + \J Q 2 y(l - y) - q\y 2 sm 2 (</>i - (j) q ). 
The lightcone momentum fractions are in terms of the lab-frame coordinates given by 



± yJQ 2 +Qt 



(22) 



(23) 



(24) 



(25) 



Having all those ingredients the lepton and hadron tensor can be calculated. In the hadron tensor the k-r and 
integrals are performed. After contracting the lepton and hadron tensor using Eq. ([7]) , the resulting expression for the 
cross section is integrated over y and expanded in powers of qr/Q except for the Gaussian in the distributions, which 
delivers the high qr suppression, and the expression for xx$ in the distribution functions. This expansion allows us 
to perform the <f> q integration analytically. After having done the <p q integration we obtain the following approximate 
expression for the symmetric cross section, 



da'- 



dqrdQd^idY 



^ 9(k 2 )Qs 



(26) 



which is accurate up to leading order in O^t/Q), furthermore we have defined 

f?{ Xi ,x 2 ) = f!( Xl )f!(x 2 ) + f!(x 2 )f!( Xl ). 

This cross section integrated over <fii is plotted as function of the three remaining variables in Fig. [2] 
For the antisymmetric cross section we find the expression, 



da 1 



dqrdQdipidY 



E 



a 2 e 2 q \S T \ 2 q T 
9M£Qs 



,-g?,/2(fc|>. 



_|_ g — 9t/ 2 (' c t)wg 



9t 



It 



2(k 2 ) s 16Q 2 (k 2 ) s 



■ cos 2<p l s 



qj, 



q T 



2{k 2 ) 



T/WG 



16Q 2 (fc2) 



■ cos 2<t> s 



T/WG 



(27) 



(28) 



G q 1T (x 1 ,x 2 ) >, 
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da/dq T dQdY [pb/GeV 2 ] da/dq T dQdY [pb/GeV 2 ] do/dq T dQdY [pb/GeV 2 ] 




(a) Q = 5 GeV and y = (b) q T = 1 GeV and Y = (c) Q = 5 GeV and q T = 1 GeV 



FIG. 2: Differential cross section for the Drell-Yan process at RHIC energy y/s = 500 GeV. 



in which we kept leading order terms in qr /Q only, which is for the </>; independent part 0{qr /Q) and for the cos 2<p l s 
dependent part 0(q^/Q 3 ). Furthermore, (fr s = (f>$ — 4>i is the angle between the spin plane and the lepton transverse 
momentum and 



F#(x u x 2 ) ee fl T q (x 1 )f^(x 2 )+f^(x 2 )f^(x 1 ), 
G\ T {x u x 2 ) = gl T (xi)gl T (x 2 ) + gl T {x 2 )g\ T {x 1 ), 



(29) 



in terms of f^(x) and g\ T (x), which are defined through the relation 



~\ K Tls 



in which 



(30) 



T/wg 



{k 2 T )Ml 
(4) + M? 



<#>. - 77^t¥t72 - (31) 



We will define the three spin asymmetries 



A TT{qT) = kJ^l^ 
J d(pida s 



A^ T (q T ) = ^ 1 , (32) 



/ r7r/4 f37r/4 f 5rr/4 f 77r/4\ / 

(Li/4 - L/l +h./4 ) d ^ da 

J 27r d(pida s 

(r /2 -r /2 +r /2 -/3 2 ; /2 ) ^ 



A TT (q T ) _ r 2rr ~ ~ ~ a ' 

J d(pida b 

which select out the <j> l s independent, cosine modulated and sine modulated terms, respectively. The latter, A^, T (q T ), 
will be zero in this case, but not in W production (cf. next section) or j-Z interference (not considered here). Both the 
A^, T (q T ) asymmetry, to which transversity contributes, and the A^ T (q T ) asymmetry receive a nonzero contribution 
from the double Sivcrs and WG effects and can be written as 



A TT {q T ) = 



! < fc T> J "J? (, 4 \Z q e 2 g F#(x 1 ,x 2 ) 



, <4>-<4> WG / 2 \ e 2 G q ^(x- i x?) 



1 2(4)woJ E q e^(x u x 2 ) j (33) 



and 



32nM*Q* \ (fc|) s E q e 2 q F?( Xl ,x 2 ) ' (k*) WG T, q e 2 q F?(x u x 2 ) 
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We note that the bound on the cos 2<p s double transverse spin asymmetry as a function of qp from transversity was 
estimated, within a collinear Collins-Soper-Sterman resummation approach to be maximally of order 5% and 
fairly flat in qp up to a few GeV at RHIC at a center of mass energy of 500 GeV. The first extraction of the quark 
transversity distribution h\ 0, however, indicates it to be about half its maximally allowed value at Q 2 ~ 2 
GeV 2 . Therefore, if this also applies to the antiquark h\, an asymmetry of 1% or less should be expected at RHIC. 

Asymmetries that are below the permille level in the entire kinematic range of interest will generally not be shown. 
They will be below the detection limit at RHIC, which will be mainly restricted by systematic errors. For the case of 
Drell-Yan this will only leave the Sivers effect contribution to the asymmetry A^ T (qp), displayed in Fig. [3] as function 
of qp, Q and Y. In the plot we also included, albeit completely negligible, the Sivers effect contribution to A^ T (qp), 
just in case the Sivers function at these values of x and Q turns out to be much larger. 

The Aj, T {qx) asymmetry reaches up to the percent level, but only for large Q 2 outside the range of interest. In the 
standard Drell-Yan range between the J/t/j and T, the asymmetry is on the permille level for the double Sivers effect 
and far below that level for the double WG effect. 

The Aj, T {q T ) asymmetry receives a contribution from the double Sivers effect at a level of 10 -6 and from the gyp 
function a contribution at a level of 10 -8 . At small Q the asymmetry is small due to the smallness of the Sivers 
function with respect to the unpolarized distribution at low values of x, whereas at higher values of Q the qp/Q 2 
suppression becomes important. One way or the other, these magnitudes are far below the detection limit at RHIC, 
even if one takes into account a possible enhancement of the effect by an order of magnitude due to the uncertainty 
in the used parametrization of the Sivers function. Therefore, the TMD effects will not spoil a determination of the 
transversity distribution if those are determined from A^ T {qp) in the lab frame instead of in the Collins-Soper frame. 
As a cross-check, to assure that TMD effects are small, one could verify that the Aj,p{qp) asymmetry is small. The 
App(qp) asymmetry is bounded by the larger Aj, T {q T ) asymmetry due to the qp/Q 2 suppression, irrespective of any 
assumptions on the Sivers function or the Worm-Gear distribution. We want to note that, considering asymmetries 
of this size, higher twist effects could become important. In case of incomplete averaging over the azimuthal angle, 
the 4>i independent asymmetry Aj.p(qp) may form a background for a determination of the (f>i dependent Aj, T (qp), 
but given its magnitude this should also not pose a problem. 




(a) Q = 5 GeV and Y = (b) q T = 1 GeV and Y = (c) Q = 5 GeV and g T = 1 GeV 



FIG. 3: Contribution to Att(qt) in the Drell-Yan process from the double Sivers effect at RHIC energy y/s = 500 GeV. 

The q'T-integrated asymmetries have also been calculated and were found to be a factor two smaller for Aj.p 
and a factor 1000 smaller for Aj. T , This agrees with the expectation that such effects are (at least) 0(M 2 /Q 2 ) 
power suppressed in this case. For completeness, we mention that the maximal qr-integrated A^ T asymmetry from 
transversity is estimated to be at the few percent level at RHIC at a center of mass energy of 200 and 500 GeV 

mm 

V. SPIN ASYMMETRIES IN W BOSON PRODUCTION 

In W-boson production one expects zero contribution from transversity within the Standard Model @, H| ■ This 
strong prediction offers the possibility to study contributions from a hypothetical W boson that appears in many 
extensions of the Standard Model. It is at least ten times heavier than the W boson and is expected to produce 
double transverse spin asymmetries at RHIC only through its mixing with the W boson, which gives the latter a 
small righthanded coupling to fermions. In Ref. [f5j this was studied quantitatively and it was shown that RHIC may 
be able to produce competitive bounds if the original design goals are met and if h\ / ff is not very much smaller than 
h\/ '/'. It was also shown that all other Standard Model background was well below the permille level, but for the 
background from the double Sivers and WG effect this was not yet quantitatively estimated. In the leptonic decay 
of a W boson the neutrino will go unobserved, which renders it impossible to determine the Collins-Soper frame and 
remove the background. Therefore, we will estimate quantitatively the sizes of these Standard Model effects to see 
whether it would jeopardize any measurement of this W — W mixing using spin asymmetries. 



9 



In Ref. [HI the TMD background was dismissed on the basis of a dimensional counting argument. If the single Sivers 
effect is a 10% effect, a double Sivers effect asymmetry in W production would be on the percent level. However, as the 
Sivers asymmetry is an azimuthal asymmetry, the lack of the knowledge of the W momentum prevents reconstruction 
of the asymmetry in W production directly. Instead, a lepton asymmetry can be measured (cf. [50]), which has a 
reduced magnitude. Naively one would expect from a dimensional analysis a large suppression of the size q T /l T , where 
qr denotes the gauge boson transverse momentum and It the lepton transverse momentum. The reason being that 
the asymmetry should vanish in the limit qT — > and the only compensating scale is It- This yields an asymmetry 
well below the permille level. A similar argument would suggest the single spin asymmetry .A/v(<7t) in W production 
arising from the Sivers effect to be qr/W suppressed, leading to an asymmetry below the percent level. However, the 
Sivers effect in A N (q T ) in W production has recently been studied theoretically [|[ and a large asymmetry (of order 
10%) was predicted. Moreover, in Ref. [T(| the lepton asymmetry A^{It) was evaluated numerically, which has a 
reduced magnitude, but still is around 3% for W + production. This is larger than expected from the dimensional 
argument and is likely because near resonance the width of the W boson becomes an important scale. The suppression 
can therefore be only as small as qr/^w- 

Similarly, a double Sivers effect contribution to Att^t) in W production is expected to be on the percent level 
and a factor q T /^w smaller for the lepton asymmetry Att{W) near resonance. When integrated over It instead, 
one can expect the asymmetry to be suppressed by a factor of q T /M^, which implies an asymmetry well below the 
permille level. Below we confirm these insights in an explicit calculation. 

Starting from Eq. , the spin symmetric and antisymmetric cross section for the production of a charged lepton 
from a W~ boson decay can be expressed as 



da s > A In 



dY T / #x tr Wg A Dp P D* wr IS° i (35) 



dl T dYid4>i 85(271-)' 

where we used qr = It +1t to write d 2 lr = cPqr- For W~ production the quark and lepton vertices are 

vf - f 2 YV L . 

The charged lepton (I) and neutrino (I) momentum 4- vectors can, in the lab frame, be expressed by 



(36) 



I = 



V2 'V2 



I = 



It -y t It y r . 

7T ! <7T < qT - lT 



(37) 



in terms of the neutrino rapidity Yj and charged lepton transverse momentum It and rapidity Y;. Having those 
ingredients, we calculate the contraction of the lepton and hadron tensor as in Eq. |[7} to get the cross section. The 
lightcone momentum fractions can be expressed in terms of It and through a power expansion in qT as 

*i = % (e Y ' + e^) ^e* cosfa - 4> q ) + O 

x 2 = ^= (e~ Y > + e-n) - ^ e -n cos (^ - cf> q ) + O ( & 

As we are working at leading twist only, we can drop the non-leading terms in this expression as well. The advantage is 
that there will not be any dependence in the distribution functions, which allows us to perform the qr integration 
in the cross section analytically. After having done the qr integration analytically, we expand in the cross section 
in parton transverse momentum up to order fc2, and p T . The integration with respect to and pt coming from 

the expression for the hadron tensor in Eq. (j4|) can now be done, which results in an expression in terms of g^,{x), 

defined in Eq. fp~6|) . and fir\x), which is likewise defined as 

f$\x) ee J d 2 k T ^f« T (x, k T ). (39) 



We find for the symmetric part of the cross section for W production 



da s _ g A \VS d KM \ 2 % 



dl T dYid(j)i 48(2tt) 2 s J 1 D' 



dY T -, (40) 
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and for the antisymmetric part 
da A 



dlxdYidfii 



g 4 M 2 \S T \ 2 \V^ M \ 2 l T f f A n B r r , „ M l 



(41) 



where 



A =150$ - Z2l%M^ - Ul^M^r 



28l^M^T 2 w 



2M^T 2 W 



44 
41% 



58$ - 9$M^ + (M 2 , + r^) - 21% (3M^ + bM^T^) cosh[y - Yj] 



26$-3M^(M, 



cosh[2(Yi - Y T )] + (24$ + 4$M^) cosh[3(F/ - Y" r )] + 2$ cosh[4(Y/ - Y T ) 



B =13(M§. + 32$M& + 16 ;2 M 4^ ^ + _ M ^ ^ + T 2 w y 



41% 



15M^ 



llM^r^) 



41% 
121% 



541% 



91%M^ 



(Afj 



2$ (9M^ + 7M^r^) cosh[l, - y r ] 



101% - Mw (Mw + r 



D =61% - 41%M ] 
and 



cosh[2(y - Y T )] + (40$ - 41%M^) cosh[3(y - Yj)] + 61% cosh[4(y - Y T )], 
(81% - 41%M? A ^ — ^ rv: ^ 1 °' 4 



i T ivi w) cosh[F/ - Yj] + 21% cosh[2(y - Yf)} 



(42) 



F° = e Y '- Y ' 



At • ( x i)fiT W ( x *) ~ Sit •(xi)9it ) (xh) 



F c 
F s 



„Y,-Y T 



„Y,-Yt 



At^i)/; 



At" WSirfci) - At ^OSit 1 (^2) 



,Yr-Y, 



„Yr-Y, 



At (x2)fiT (1 \xi) - gfrH^girH^) 



ftr ( x 2)9tr fa) - At (si)Sit faa) 



(43) 



With the use of the expressions for the cross section in Eqs. (|4"0"|) and (|4"Tj) . the spin asymmetries, as defined in Eq. 
(JH3> , can be written as 



4° 



TT 



(It) 



\S T \ 2 M 2 JdY I± %F° 



A T 'Jji (i T ) 



\S T \*M$f dY T &F™ 



(44) 



nl% 



The results are easily modified for W + production by substituting u — > d, d — > u in Eq. (1431) , substituting It — > 
Itt4>i — > 0r m a ^ expressions and integrating over Y\ instead of Yj in the cross sections and asymmetries. As a cross- 
check of the approximation method employed here, we calculated the single spin asymmetry An in W production and 
found reasonable agreement with the results in Refs. 0, [l(| taking into account that different functions were used. 
The cross sections for W ± production are plotted in Fig. 01 We only show the asymmetries in W + production 



dajdlrdYi [pb/GeV] 

- W+ 

- w- 



da/dl T dY t [pb/GeV] 




It [GeV] 



/ 4 




w 


/ 3 




w 


/ 2 
/ 1 






-1.5 -1.0 -0.5 


0.5 1.0 


.5 


(b) l T 


= 40 GeV 





Y, 



FIG. 4: Differential cross section for W boson production at RHIC energy y/s = 500 GeV. 



in Fig. [5J because they are largest. The maximal asymmetry is near resonance and reaches up to 0.15%, which is 
already below the detection limit at RHIC. However, for a bound on a possible W—W mixing it is not the differential 
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asymmetry that is relevant, but the asymmetry in the integrated cross section. In those asymmetries the contribution 
at It < Mw/2 largely cancels the contribution at It > Myy/2, resulting in very small asymmetries. We find the 
integrated asymmetry in W" production around 10 -7 and in W + production around 10~ 6 , far below detection limits 
at RHIC. This confirms the expectation expressed in Ref. 15] that the background from TMDs, is indeed negligible. 



0.0015 - 
0.0010 1 
0.0005 1 



-0.0005 - 
-0.0010 1 
-0.0015 1 




■ - A% T x 3 (WG) - - A° TT x 3 (WG) 

(a) Y t = (b) l T = 40 GeV 

FIG. 5: Contributions to Att(It) in W + boson production from the double Sivers and Worm Gear effects at RHIC energy 
Vs = 500 GeV. 



VI. SUMMARY AND CONCLUSIONS 



We calculated the transverse momentum dependent double transverse spin asymmetries in the laboratory frame 
for Drell-Yan and W production arising from the Sivers effect and from the Worm Gear distribution function gyr 
within transverse momentum dependent factorization. Those asymmetries were previously calculated only in the 
Collins-Soper frame where they are independent of the lepton azimuthal angle. The advantage being that one can, 
in that frame, easily distinguish them from the asymmetry coming from transversity. In the lab frame a residual 
dependence on the lepton azimuthal angle of the TMD asymmetries survives and enters the double transverse spin 
asymmetry in exactly the same way as transversity does. This is in contrast to a collinear factorization approach 
where the effects from TMDs are absent to begin with. Therefore, a nonzero cos 2<p l s asymmetry Att(qt) in Drell-Yan 
in the lab frame is a priori not a sufficient indication of a nonzero transversity distribution. However, from what is 
known about the magnitudes of the Sivers and Worm Gear functions, our conclusion is that the TMD background is 
below the permille level. Therefore, a percent level asymmetry can be viewed as coming from transversity. That is 
an important conclusion for the RHIC spin program. Transversity distributions can thus safely be determined from 
the transverse momentum dependent double spin asymmetry in the lab frame, like for the (^-integrated asymmetry, 
assuming of course the antiquark transversity distributions are sufficiently large. As a cross-check of the smallness 
of the TMD background, one can verify that the angular independent A TT (qx) asymmetry that arises only from the 
mentioned TMD effects, is indeed much smaller. 

In the leptonic decay of a W boson the neutrino will go unobserved, which renders it impossible to determine the 
Collins-Soper frame. In that frame a correlation between the lepton azimuthal angle and the proton transverse spin 
direction can solely be caused by a non-zero righthanded coupling of the W boson in combination with a non-zero 
transversity distribution, which makes it a very suitable process for the determination of a possible W — W mixing 
as discussed in [lSf . In the lab frame, however, there might again be a residual asymmetry coming from the double 
Sivers or WG effects. They can lead to a nonzero result in W production, which could be mistaken for physics beyond 
the Standard Model or simply spoil the opportunity to bound a possible W — W mixing. We obtained numerical 
estimates for the sizes of the asymmetries at RHIC and found that they are far below the detection limits. This 
means that the background from the double Sivers and Worm Gear effects is negligible and does not hamper the 
investigation of the complex mixing of W-W bosons as discussed in [l5| . 
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